Determination of ground state properties in quantum spin systems 
by single qubit unitary operations and entanglement excitation energies 
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We introduce a method for analyzing ground state properties of quantum many body systems, 
based on the characterization of separabihty and entanglement by single subsystem unitary opera- 
tions. We apply the method to the study of the ground state structure of several interacting spin-1/2 
models, described by Hamiltonians with different degrees of symmetry. We show that the approach 
based on single qubit unitary operations allows to introduce "entanglement excitation energies", 
a set of observables that can characterize ground state properties, including the quantification of 
single-site entanglement and the determination of quantum critical points. The formalism allows 
to identify the existence and location of factorization points, and a purely quantum "transition 
of entanglement" that occurs at the approach of factorization. This kind of quantum transition is 
characterized by a diverging ratio of excitation energies associated to single-qubit unitary operations. 

PACS numbers: 03.65.Ca, 03.67.Mn, 73.43.Nq, 75.10.Jm 



I. INTRODUCTION 

With the advent of quantum information theory, en- 
tanglement has been recognized as a fundamental phys- 
ical resource, on equal footing with other fundamental 
resources such as energy and entropy. The relationships 
between these resources are being actively investigated, 
and it has been suggested that the entanglement prop- 
erties in many-body systems may be characterized by 
suitably defined energy observables [l|, . Establishing 
direct connections between energy observables and en- 
tanglement is of interest because the former could lead 
to a deeper conceptual understanding of the latter and be 
exploited for its experimental production and manipula- 
tion in complex quantum systems. On the other hand, 
the study of the role played by entanglement in quantum 
phase transitions |3, 01 and its relations with ground 
state (GS) properties has sparked a rapidly growing field 
of research Q. As entanglement plays a fundamental 
role in quantum information Q , the question then arises 
whether it is possible to determine structural aspects 
of complex quantum systems using concepts and tech- 
niques of quantum information theory. Relevant open 
challenges include the understanding of the relation be- 
tween GS entanglement, physical observables, and criti- 
cality, as well as the determination of factorization points 
at which quantum ground states become separable [1, Q . 

Recently, a formalism of single subsystem unitary op- 
erations has been introduced [HI : These operations 
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are defined as unitary and traceless local transformations 
with non degenerate spectrum (as well as Hermitian, in 
the case of qubits) that act on a single subsystem, leav- 
ing unaffected the remaining ones. The corresponding 
approach provides a tool for studying pure state entan- 
glement for systems of qubits (spin 1/2) or qutrits (spin 
1) , an d for Gaussian states of continuous variable sys- 
tems |ll| . For such systems, a necessary and sufficient 
condition for pure state separability is the existence of 
a unique transformation (termed invariant) that leaves 
the state unchanged. On the other hand, if a pure state 
is entangled, the minimum Euclidean distance in Hilbert 
space from the state and its image under the action of 
single subsystem unitary operations singles out a unique, 
extremal operation that is directly related to the entan- 
glement properties of the state. Namely, the squared 
minimum distance coincides exactly with the linear en- 
tropy, and is thus monotone in the entropy of entangle- 
ment of the state. As a consequence of these results, 
separability points are determined and singled out by im- 
posing equality of the extremal and invariant operations. 

The approach based on single subsystem unitary op- 
erations allows to gain informations on the global nature 
of pure states of a composite system by looking at how 
states transform under local operations that, by defini- 
tion, do not change the content of entanglement (local 
unitaries). Moreover, the method finds a natural oper- 
ational interpretation, as looking at the response to a 
given action is a basic tool in the investigation of physi- 
cal properties. 

The structure of the paper is as follows: In Sec. |TT] 
we introduce Single Qubit Unitary Operations (SQUOs), 
and show that each SQUO singles out a direction in the 
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three-dimensional single-spin space; we then define the 
extremal SQUO, as the one associated to the distance 
of the ground state (GS) from the nearest state in the 
set of states obtained from the GS under the action of 
SQUOs. In Sec. IIIIl we define the excitation energy rel- 
ative to each SQUO, and demonstrate that, under rather 
general conditions, the excitation energy associated with 
the extremal SQUO vanishes if and only if the GS is fac- 
torized: We therefore name it entanglement excitation 
energy (EXE). Sec. [IV] is devoted to the analysis of sev- 
eral 5=1/2 one-dimensional spin models, with attention 
focused on the dependence of the EXE on the Hamilto- 
nian parameters. In Sec. |V] we present and discuss our 
results: we find that the EXE follows monotonically the 
behavior of the GS single-site entanglement (i.e. the von 
Neumann block entropy between a single spin and the 
rest of the system) and exhibits signatures of the rele- 
vant GS properties, including the existence and location 
of factorization and quantum critical points. Finally, we 
introduce the orthogonal SQUOs as those defined by two 
directions orthogonal to each other, and to the direction 
selected by the extremal SQUO: We find that the cor- 
responding excitation energies coincide when the EXE 
vanishes, in such a way that the ratio between their dif- 
ference and the EXE itself diverges at the approach of a 
factorization point in a large class of models: such a di- 
vergence may therefore be exploited in order to detect the 
entanglement transition associated with the occurrence of 
a fully separable ground state p^, Ts^l in strongly corre- 
lated quantum systems. In Sec. IVII we draw our conclu- 
sions and discuss some possible future lines of research. 

II. SINGLE QUBIT UNITARY OPERATIONS 

Consider a iV-qubit system. A SQUO is defined 
as the following unitary transformation 

Uk = (g) 1, «) Ofc , (1) 

where the operators 1^ are the identities on the N ~ 1 
spins different from spin fc, and the operator Ok is a 
Hermitian, unitary, and traccless 2x2 transformation, 
acting on spin k. In the standard basis {| j), | J,)} , Ok 
is parametrized as 

where 9 and (p vary in the ranges [— 7r/2, tt/2) and [0, 27r), 
respectively. The above expression may be written as 
Ok = u ■ (Tk, where CTk = i'^k^'^k^'^k) Pauh 
matrices, and the unitary vector u = {ux,Uy,Uz) = 
(sin 9 cos (fi, sin 9 sin (p, cos 9) defines a direction in the 
three-dimensional single-spin space: each direction in 
such a space defines a SQUO and, viceversa, to each 
SQUO there corresponds one such direction. The ac- 
tion of a SQUO transforms any pure state j'l') defined in 



the 2^ dimensional Hilbert space in the state Uk\'^)- If 
C/fcl^') = I^E") then Uk is the invariant SQUO. In general, 
however, the transformed state differs from I^P): This 
difference can be quantified by introducing the standard 
Hilbert-Schmidt distance: 



d(C/,;|*)) EE v/l-|(*|C/fe|*)P = 

= ^l-mu-ak\^f)\^^diu;\M>)) . (3) 

Minimizing the distance over the entire set of SQUOs, i.e. 
evaluating minje^^} d{u; l^*)), yields that the distance of a 
pure state from the nearest transformed state is obtained 
for [13: 

,~=arctanfMV 

9 = arctan ^Ml^^^l^M^^ , (4) 

where denotes the spin expectation values {(j^)/2 
{a = X, y, z) on the state j^*). The values ^ fix the direc- 
tion u which determines the extremal SQUO, C/|^*''. The 
distance d{u\ Y^)) of a state from the nearest transformed 
state, corresponding to the extremal SQUO U^^*^ , can be 
evaluated explicitly, and its square reads 

d'{u; m) = 1 - [(A/f )2 + (M|)2 + (M|)2] . (5) 

This is exactly the linear entropy Sl of state j^*): Sl = 
2(1 — Tr[p|]), where pk is the single-site reduced density 
matrix after tracing out the — 1 remaining spins. It 
is a well known fact that in the case of qubits the linear 
entropy coincides with the tangle r: Sl = t = 4Detpfe 
[l3,[l3|- The pure state entanglement, measured by the 
von Neumann entropy £ [l^, is a single- valued, mono- 
tonic function of either of these two coinciding quanti- 
ties (either the hnear entropy or the tangle): = 
~xln2X — (1 — x)ln2(l — x) where a; = (1 + Vl — t)/'^- 
Therefore, due to Eq. ([5]), extremal SQUOs and the as- 
sociated Euclidean distances determine pure state entan- 
glement, whose entropic quantification is recast in direct 
geometric terms. 



III. ENTANGLEMENT EXCITATION 
ENERGIES 

The full consequences of the above results can be ex- 
ploited for the characterization of GS properties. In the 
following, we show that t/^^*'' can be used to construct 
an energy observable which witnesses GS separability and 
quantifies GS entanglement. 

Let \G) be the GS of a system of interacting spins with 
Hamiltonian 7i, and define the excitation energies, 

AE{Uk) EE {G\ulHUk\G) - {G\n\G) ; (6) 
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the dependence of AE{Uk) on the parameters 9 and (p, 
i.e. on the direction u selected by the SQUO Uk, is made 
evident by using Eqs. ([T]) and One obtains that 

AE{Uk) = ulAE{al) + ulAE{al) + ulAE{al) + 

where Eq^/j = ep^a = (G|CTj?Ha-f + a'^Ha'j^\G). 

Amongst these excitation energies, which are non- 
negative by definition, the one corresponding to the ex- 
tremal SQUO, A£'([/^^*''), has a prominent role in the 
analysis of entanglement properties of the GS, due to 
the intimate relation between its vanishing and the oc- 
currence of full separability. We will therefore name it 
entanglement excitation energy (EXE). 

In fact, if in the GS the spin k is not entangled with 
the rest of the system, then Uj:^*^\G) = \G), and the 
EXE consequently vanishes. On the other hand, if the 
spin k is entangled with the rest of the system, one nec- 
essarily has Uji^^^lG) ^ \G) and the EXE may vanish 
only if U^^*^\G) is again a state of minimum energy, i.e. 
if (G|[H,J7^^*'']|G) = 0. Therefore, barring pathological 
or trivial cases in which the Hamiltonian does commute 
with each SQUO, the vanishing of the EXE occurs only 
in the presence of a factorized GS, proving the statement. 

It is worth noticing that the vanishing of the EXE is a 
necessary, but not sufficient, condition for the separabil- 
ity of excited states. In fact, considering an entangled ex- 
cited state \tp), the action of the corresponding extremal 
SQUO transforms it in a different state that is not, in 
general, an eigenstate of the Hamiltonian and can have 
a non vanishing projection on any possible eigenstate of 
the system. The possibility to populate eigenstates with 
energy lower than that associated to \ip) makes then it 
possible to obtain a vanishing EXE even in the presence 
of an extremal SQUO. 

IV. INTERACTING SPIN SYSTEMS 

Equipped with these results, we move on to apply them 
to the determination of the GS properties of interacting 
spin systems. We consider antiferromagnetic anisotropic 
Heisenberg-like XYZ model Hamiltonians: 

= ^ SfSf:^i + AySfSf^-^ + A^S-S^^-^ - hS- , (8) 

i 

where Sf — cr°'/'2, the index i runs on the sites of a 
one-dimensional lattice, and periodic boundary condi- 
tions are assumed. In Eq. ([S]) /i is the reduced (dimen- 
sionless) magnetic field, and Aj,, are the anisotropy 
parameters, < Ay^z < 1- This choice is motivated 
by its great generality: according to the different val- 
ues of the anisotropics, Eq. ^ comprises many differ- 
ent models belonging to different classes of universal- 
ity, including the Ising and the Hcisenberg ones. For 
the models described by Eq. ([S]) genuine quantum phase 
transitions, with or without order parameter, occur at 



zero temperature at critical values he that depend on the 
anisotropics. If Ay = 1, the symmetry forces a transi- 
tion without order parameter, and the system remains in 
a critical regime for all h < h^. In this case if A^ ^ 1 
the system is in the XX Z symmetry class and goes over 
to the Hcisenberg symmetry class when A^ = 1. For 
< A J, < 1, the system undergoes a second order phase 
transition, acquiring a non- vanishing staggered magneti- 
zation along the x direction for h < he- Antiferromag- 
netic order divides the lattice in two sublattices, each 
characterized by opposite value of the magnetization. 
Again, changing the value of A^ lets the system move 
from the Ising symmetry class (A^ = 0) to the XYZ 
symmetry class (A^ ^ 0). Concerning the magnetization 
along the directions y and z, one has that, regardless of 
the anisotropy, the magnetization My is always vanish- 
ing, while AIz vanishes at and only at h = 0. Importantly, 
in the models described by Eq. ([5]) the external field not 
only drives the system through a quantum phase transi- 
tion, but induces as well the factorization of the GS, for 
h ^ hf = y/{l + Az){Ay + Az) This phenomenon 
has been recently investigated in relation to the analysis 
of entanglement properties in spin systems [§]. For the 
ground states of the models described by Eq. ([8|) , the van- 
ishing of the EXE is a necessary and sufficient condition 
for GS separability. 

Using Eqs. ([4]), and considering that it is always — 
0, and hence = 0, on gets, by Eq. ([7]), the exact ex- 
pression of the EXE that reads 

AE{u) = A£;(C/'=^*'') = 

_ 9xxMl - MxMzgzx ^ ^yUyy 

, Az [gzzMl + MxMzgzx) ] 



where r = 5^ is the GS tangle (linear entropy), gap = 
(GIS'f S'f|_j|G) are the nearest-neighbor correlation func- 
tions. Ma — (Gl^f |G) are the expectation values of the 
spin operators, and any dependence on the site index is 
neglected due to translational invariancc. 

V. RESULTS 

We first analyze the EXE, AE{u), and compare it with 
the von Neumann entropy £. The latter measures the 
bipartite single-site entanglement between one selected 
spin and the rest of the chain , and it provides an up- 
per bound to all bipartite block entanglements. There- 
fore, its vanishing guarantees the full separability of the 
GS. Exploiting the conditions for the vanishing of Sl and 
AE{u) allows to determine unambiguously the value and 
location of the factorization point. 

In Figs.[T]and[2]wc plot AE{u) and S as functions of h. 
Fig. [1] shows exact analytical results [13, [3 the XY 
model, and Quantum Monte Carlo (QMC) data for the 
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FIG. 1: (color online) Entanglement excitation energy AE{u) 
(Black solid line) and single-site entanglement £ (Red dashed 
line) as functions of the reduced field h, for a): Ay = 0.4, 
= (XY), and b): Ay = 0.25, = 1 (XYZ). AU 
quantities are dimensionless. 

XYZ model. Fig. [2] shows exact analytical results for 
the XX model; and QMC data for the Heiscnberg model. 
Numerical data have been obtained via stochastic-series- 
expansions quantum Monte Carlo simulations [l9j . 
based on a modified directed-loop algorithm. The EXE 
and the entropy of entanglement evidently exhibit the 
same behavior, vanishing simultaneously at the factor- 
ization points. At critical points, the derivatives of both 
AE(u) and £, with respect to h, diverge, signaling the 
onset of a quantum phase transition. This divergence 
occurs both for = 1 and < 1, at variance with 
the concurrence, i.e. a measures of the entanglement be- 
tween two spins of the lattice [2^, that exhibit different 
behaviors at the critical points for models with different 
symmetries In this sense, we suggest that the EXE, 
contrary to other observable estimators of entanglement, 
can be considered a universal indicator for the onset of 
quantum phase transitions in spin systems. 

The relevance of SQUOs and EXEs in determin- 
ing entanglement and factorization properties, and in 
general qualitative changes in the GS can be under- 
stood in physical terms and related to "entanglement 
transitions" [H, [3 ■ If a system is in a pure classical 
state the orientation of any spin is well defined. Select- 
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FIG. 2: (color online) Entanglement excitation energy AE{u) 
(Black solid line) and single-site entanglement £ (Red dashed 
line) as functions of the reduced field h, for a): Ay = 1, 
A^ = (XXZ), and b): Ay = 1, A^ = 1 (Heisenberg). All 
quantities are dimensionless. 



ing a spin and performing a rotation about its orientation 
leaves the state of the system and its energy unchanged. 
Viceversa, if the system is in an entangled quantum state 
the single-spin orientation is not defined and any "rota- 
tion" , i.e. any SQUO, will change the state of the system 
and will yield a change in energy. Therefore, the increase 
in energy is connected to "how much the GS is entan- 
gled" . Constructing the excitation energy associated to 
the extremal SQUO, i.e. the EXE, formalizes the argu- 
ment. 

Based on the above comparison with the classical case, 
let us introduce two directions, u'j_ and m" , orthogonal 
to each other and to the direction u associated to the 
extremal SQUO f/'^^*''. Being u = (sin6i, 0, cosS), we 
choose u'^ — (cos 6', 0, — sin 6*) and it" = (0,1,0). These 
directions define the orthogonal SQUOs, via Eq. and 
the corresponding excitation energies AE'j_ and Ai?" , via 
Eq. (I7|). In Fig. Owe compare the field dependence of the 
EXE with that of the above defined excitation energies, 
for the XY and XYZ models: We find that AE'j^ and 
AE'^ coincide at hf, i.e. when the GS is factorized, in 
full analogy with the classical case. In fact if the system 
is in a pure classical state any rotation of 7r/2 around a 
direction orthogonal to the spin orientation causes an in- 
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FIG. 3: (color online) Entanglement excitation energy AE{u) 
(Black solid line), excitation energy AE'j_ (Red dashed line), 
and excitation energy AS" (Blue dotted line) as functions 
of h, for a): Ay = 0.4, = [XY), and b): Ay = 0.25, 
A^ = 1 (XYZ). All quantities are dimensionless. 



crease in energy by one (in units of the reduced field). For 
h hf, AE'j^ and Ai?" behave differently, a fact which 
has no classical analogue, so that the deviation from clas- 
sicality can be quantified by the difference AE'j_ — AE'^ . 

We can then compare such difference to the amount of 
entanglement, as measured by the EXE, by defining the 
entanglement energy ratio (EER) 



Re — 



AE'^ - AE'l 
AE{u) 



(10) 



In Fig. m we show the behavior of the EER Re, as a func- 
tion of /i, for two models belonging, respectively, to the 
XY and to the XY Z universality class. In both cases, 
Re diverges at the factorization point hf, and its first 
derivative diverges at the critical point he- The abrupt 
change in the GS properties at the factorization point, as 
measured by the divergence of the EER, signals a quali- 
tative change of purely quantum nature, a "transition of 
entanglement" that occurs when hf is approached. We 
remark that all the conventional on-site expectations and 
n-point correlation functions remain finite and analytic 
at a factorization point, confirming that the divergence 
of the EER cannot be associated to any kind of conven- 
tional quantum phase transition. On the other hand, the 
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FIG. 4: (color online) Entanglement energy ratio Re as a 
function of h, for a): Ay = 0.4, A^ = {XY), and h): Ay = 
0.25, Az — 1 (XYZ). All quantities are dimensionless. 



EER remains fixed at the constant value Re = —1 if we 
consider models of interacting spins with = 0: The 
fact that in this case the EER fails to identify the tran- 
sition from entanglement to factorizability is intriguing 
and deserves further studies. 



VI. CONCLUSIONS AND OUTLOOK 

In summary, we have defined entanglement excita- 
tion energies (EXEs) associated to the extremal single 
qubit unitary operations (SQUOs). We have showed that 
EXEs are useful tools for the study of various ground 
state properties for interacting spin systems, including 
the determination of factorizability, the quantification of 
single-site entanglement, and the identification of quan- 
tum critical points. We have discussed how SQUOs and 
EXEs determine unambiguously the presence or the ab- 
sence of a separable GS, and we have introduced an 
entanglement energy ratio (EER) of excitation energies 
that diverges (for Aj, ^ 0) at the approach of a factor- 
ization point, thus defining an entanglement-separability 
transition of purely quantum origin. 

Considering possible future developments along this 
line of research, we would like to observe that when 
moving away from the factorization point, the physics 
of quantum ground states can in principle be recovered 
by expanding the single-site entanglement in powers of 
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the EXE, as the latter naturally defines a small expan- 
sion parameter of the theory around a factorization point. 
The analytic generalization of the formalism to general 
spin models is a subject of current investigation. The 
results of such an analysis could allow to determine ex- 
act solutions of non exactly solvable models at factoriza- 
tion points, as well as the approximate description of GS 
properties out of the factorization point by controlled ex- 
pansions in powers of the EXEs. The necessary and suf- 
ficient conditions for factorizability obtained within the 
formalism of SQUOs might be extended to the case of 
models involving interacting systems of arbitrary local 
dimension. In particular, it would be very interesting 
to consider interacting systems of very high spins and 
study them in a continuous-variable representation, for 
which single-subsystem unitary operations can be read- 
ily expressed in terms of single-mode unitary transforma- 
tions. This line of investigation could then be extended 
as well to include models of interacting harmonic and 
anharmonic chains and lattices. 

Conceptually, the present work discloses an intimate 
connection between two different universal physical re- 
sources, energy and entanglement, and might have prac- 
tical consequences for the experimental production and 
manipulation of entanglement, and information transfer 
in real systems of interacting qubits. The single-qubit 
excitation energy provides a well-defined form of macro- 
scopic observable for the detection and determination 
of single-site entanglement in many-body systems, along 
lines close in spirit to pair-entanglement detection in sys- 



tems of interacting magnetic dipoles by measurements of 
heat capacities and magnetic susceptibilities [2l|. As a 
working example, the method has been applied and illus- 
trated in quantum spin- 1/2 models, that are of particular 
relevance both for quantum information and condensed 
matter physics. However, in principle it can be applied as 
well to more general instances [10, UM i such as systems 
of interacting qutrits, Hubbard models, and harmonic 
lattices of continuous variables. The choice of nearest- 
neighbor couplings in the case of interacting qubits de- 
serves a comment. In fact, the method is in no way lim- 
ited by the choice of the interaction. Generalizations to 
spin models with interactions of arbitrary range and lat- 
tices of different topologies are possible and, as suggested 
above, can be of particular relevance in establishing the 
existence of factorization points and in constructing con- 
sistent descriptions of GS quantum physics by systematic 
expansions in powers of the EXE around the factorized 
solutions. 
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